Inter-site Coulomb interaction and Heisenberg exchange 
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—A of the d-orbitals is taken to be negative, and in ad- 
dition to the d-site Coulomb repulsion U we include an 
additional Coulomb repulsion V between nearest neigh- 
bor d and ligand orbitals. 

At 'half-filling' and in the limit V— 0, [/— >oo we have 
precisely one electron per lattice site for the Hubbard 
model, and one electron (or hole) per d-site in the charge 
transfer system (see the states labeled '0' in Figure [l]). 
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Based on exact diagonalization results for small clusters we discuss the effect of inter-site Coulomb 
repulsion in Mott-Hubbard or charge transfers insulators. Whereas the exchange constant J for 
direct exchange is enhanced by inter-site Coulomb interaction, that for superexchange is suppressed. 
The enhancement of J in the single-band models holds up to the critical value for the charge density 
wave (CDW) instability, thus opening the way for large values of J. Single-band Hubbard models 
with sufficiently strong inter-site repulsion to be near a CDW instability thus may provide physical 
realizations of t— J-like models with the 'unphysical' parameter ratio J/i«l. 
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A frequently used approach to discuss the strong- 
coupling limit of Hubbard-type models is the pertur- 
bation expansion in t/U EtB, which leads to an ef- 
fective Hamiltonian operating in the subspace of states 
without double occupancies, the well-known t— J model. 
By virtue of its perturbational derivation, the relation 
J = 4t 2 /U holds, so that it seems that a necessary con- 
dition for the validity of <— J-like models is the relation 
J/t = 4(t/U) <C 1. In this work we show that this is not 
necessarily the case; rather, Hubbard models with an ad- 
ditional inter-site repulsion which is sufficiently strong 
to bring the system close to a charge-density- wave insta- 
bility may provide 'physical' realizations of t—J models 
with J/t w 1. We also show that the situation is reversed 
in a charge-transfer insulator, where inter-site Coulomb 
repulsion on the contrary suppresses the superexchange. 
To begin with, it is advantageous to develop a simple 
picture of a correlated insulator. We restrict ourselves 
to ID systems in all what follows, the main reason be- 
ing that exact diagonalization for reasonably large 2D 
clusters of the single or even two-band Hubbard model 
is not feasible. We expect, however, that the physical 
considerations outlined below retain their validity in any 
dimension. We consider the single band Hubbard model 

H = -ty^( cl^Ct+La+H.c )+U ni^m,i+V nim + i 

i,(T i i 

where ni. CT =cj CT Ci. CT , »j=?ii,t+?ij,.p The model already in- 
cludes an additional repulsion between electrons on near- 
est neighbors ~V. We also consider the charge tranfcr 
system 

H=-t J2 ( d L c j,° + H -c) + U Y1 n ».T«U 

which consists of strongly correlated 'cf-orbitals sepa- 
rated by uncorrelated ligands (N(i) denotes the two lig- 
and sites neighboring the d-site i). The on-site energy 
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FIG. f. Charge fluctuations in the single band Hubbard 
model (a) and the charge transfer system (b). The dark cir- 
cles in part (b) indicate the 'd' sites, the light ones the ligands. 
The lower part of the figure gives the energies of the states cre- 
ated by the charge fluctuation with (dashed line) and without 
(full line) a nearest neighbor repulsion of strength V taking 
the hopping integral to be 0. 



Switching on the inter-site kinetic energy then will lead 
to charge fluctuations: in the Hubbard model, an electron 
can jump to a nearest neighbor, thus creating a 'hole' and 
a 'double occupancy' on nearest neighbors. The energy 
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thereby is increased by the large amount U. In subse- 
quent steps hole and the double occupancy can separate 
even further, but all states generated in this way have 
the same large Coulomb energy U (see Figure |l|) . For the 
case of the charge-transfer system, similar considerations 
hold, with the energy required to 'drag apart' hole and 
double occupancy being min(A,U). By analogy with 
the problem of a particle moving in a ID attractive 5- 
potential one may therefore expect that the probability 
to find a hole and a double occupancy at a distance of 
n lattice spacings decreases rapidly with n. The canon- 
ical transformation to a Heisenberg or t—J Hamiltonian 
with only nearest neighbor exchange then obviously is 
equivalent to assuming that states where hole and dou- 
ble occupancy are more distant than nearest neighbors 
(nearest neighbor d-sites for the case of the charge trans- 
fer system) can be completely neglected. 
Let us now consider the case V > 0. For the single-band 
Hubbard model the energy necessary to create a charge 
fluctuation on nearest neighbors is lowered to U — V (see 
Figure |l|) . The energy to drag apart hole and double oc- 
cupancy, however, remains unchanged, so that the sepa- 
ration between the two Hubbard bands (which are associ- 
ated with 'free' propagation of hole or double occupance) 
stays unaffected ||f|]. Invoking again the analogy with 
a particle moving in a ID attractive potential one might 
instead expect the formation of a 'bound state' between a 
hole and a double occupancy on nearest neighbors, which 
manifests itself as an excitonic excitation at an energy 
well below the Hubbard gap in the optical conductiv- 
ity Q. The V-term. therefore does not reduce the Hub- 
bard gap, and in particular does not 'screen' U to turn 
the correlated insulator into a weakly correlated metal. 
The question remains, however, whether a transforma- 
tion to a t—J model with only nearest neighbor exchange 
still makes sense in the presence of the y-term. Obvi- 
ously, if there is an appreciable probability to have the 
hole-double occupancy pair on more distant than near- 
est neighbor sites, longer range spin correlations will be- 
come important in the effective low energy Hamiltonian. 
To address this question we use exact diagonalization of 
small clusters; we note that for the short-range/high en- 
ergy processes which mediate the Heisenberg exchange 
finite size effects probably play only a minor role. We 
introduce n ijd = ni^Ujj and m t h = (1 - n 4 ,|)(l - n^i), 
i.e., the density operators for doubly occupied and empty 
sites and study the 'empty-double correlation function' 

9{R) = J2(n t jin l+R . d ) (1) 

i 

(where N denotes the chain length) as a function of the 
nearest neighbor repulsion V. This is shown in Figure |^. 
Obviously g(R) is peaked at R = 1, i.e., independently 
of the value of V 'hole' and 'double occupancy' strongly 
prefer to be on nearest neighbors. For R > 1 and not too 



large V the empty-double correlation immediately lev- 
els off to a practically i?-independent value goo] this can 
be understood as the contribution from holes and double 
occupancies belonging to different 'pairs': since the 
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FIG. 2. (a) 'Empty- double correlation function' g(R) in the 
ground state of a 14-site ID Hubbard ring for different values 
of V. The other parameters are (7=10, t=l. 
(b) Asymptotic value g^ vs. P2, the probability to have two 
double occupancies in the ground state. 

positions of the pairs are uncorrelated, this component 
has no more i?-dependence. To check this interpretation, 
Fig. ||d shows goo plotted versus P2, the probability to 
have precisely two double occupancies in the system (for 
our relatively short chain and large value of U the proba- 
bility to have n double occupancies decays exponentially 
with n, so that we may safely neglect corrections due to 
states with 3 or more double occupancies). To good ap- 
proximation, N ■ goo equals the probability to have more 
than one double occupancy in the chain, consistent with 
our interpretation. Assuming that a Heisenberg antifer- 
romagnet with exchange constant J = 4t 2 /(U — V) de- 
scribes the system well, we can moreover give a simple es- 
timate of g(R = 1): the number of doubly occupied sites 
is (N d ) = du-yEo, whereto = - N ■ (At 2 / (U — V))ln 2 is 
the ground state energy of the Heisenberg antiferromag- 
net. Taking into account that on the average 50% of the 
double occupancies are to the right of 'their' hole we find 
Ng(R = l) = (l/2)-(N d )=2N(t/(U-V)) 2 In 2. With the 
parameters used in Figure we find e.g. Ng(R = 1) = 
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0.1664 for V = 0, Ng(R = 1) = 0.6654 for V = 5 which 
values obviously provide reasonable estimates. 
For large values of V ~ C//2, i.e. the critical value for 
the onset of the CDW instability in ID [|| we see in ad- 
dition a modulation of g{R) at longer distances with a 
period of two lattice spacings. This is obviously is a pre- 
cursor of the CDW ordering which occurs for larger V. 
We may thus expect that V = U/2 represents a limiting 
case, where the mapping to a nearest neighbor-exchange 
model starts to become inaccurate. 

The empty-double correlation function g(R) thus demon- 
strates that the main effect of V is the enhancement of 
nearest neighbor charge fluctuations. This suggests that 
simple nearest neighbor Heisenberg exchange remains a 
good approximation also for the case V>0, the main dif- 
ference as compared to the case V — being the enhance- 
ment of charge fluctuations and hence the Heisenberg ex- 
change J. To check this we consider the dynamical spin 
correlation function S(q,u), shown in Figure ^a for dif- 
ferent values of V. 
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cally identical, as demonstrated in Figure ||b. Assuming 
the relation J = At 2 jU to hold for V=0, we can use the 
rescaling factors to extract the effective Heisenberg ex- 
change J e ff(V), shown in Figure ||c. Obviously values 
of J e ff/t, as large as 0.8, i.e. far outside the 'pertur- 
bational' range are possible. Figure ||c also shows the 
simple estimate J e ff(V) — At 2 /(U — V), which provides 
a reasonable approximation for the numerical values. Let 
us stress that for small values of V this functional form 
of J e ff(V) is rather trivial; what is not, however, is the 
validity of this relation up to the value V=U/2. 
We briefly turn to the charge transfer system. From Fig- 
ure [l] we may expect that the V-teim effectively increases 
the charge transfer energy A to A + V (it should be 
stressed, that this refers to 'bare' or 'model' parameters; 
in the case of the charge transfer insulator, the energy 
required to drag apart hole and double occupancy 
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FIG. 3. (a) Spin correlation functions of a 10-site ID Hub- 
bard ring for different V (U=W, t=l). 

(b) Same as (a) but with the frequencies rescaled so as to 
make the energies of the dominant peak at 27r/5 identical for 
all spectra. 

(c) Effective Heisenberg exchange, Jeff(V), calculated from 
the rescaling factors used in (b), the dashed line shows 
J eff (V)=4t 2 /(U-V). 

With increasing V, all peaks move to higher energy, the 
overall shape of the spectra being unaffected. Rescaling 
the frequencies the different S(q,u>) can be made practi- 



FIG. 4. (a) Spin correlation functions of a 6-unit cell ID 
Charge transfer ring for different V (A=3, (7=8, t=l). 

(b) Same as (a) but with the frequencies rescaled so as to 
make the energies of the dominant peak at 27r/3 identical for 
all spectra. 

(c) Effective Heisenberg exchange, J e ff(V), calculated from 
the rescaling factors used in (b). The dashed line has been 
calculated from 4 th order perturbation theory. 

becomes min(A + 2V, U) in the presence of the near- 
est neighbor repulsion, so that the experimentally mea- 
sured charge transfer gap would be increased to A + 2V) . 
Contrary to the case of direct exchange, this suggests a 
suppression of the exchange constant J. This is demon- 
strated in Figure |], which shows the dynamical spin cor- 
relation function for the charge transfer system with dif- 
ferent values of V. From the rescaling of the SCF (see 
Figure |^b) we can again extract the effective exchange 
constant, which is shown in Figure 0c. The figure also 
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shows a simple estimate based on 4 th order perturbation 
theory, where J={2t pd /(A + V)) 2 ■ (1/A + l/U), which 
however provides only a rough estimate. 
We return to the single band model, and consider the 
doped case. In this case, the canonical transformation 
gives the strong-coupling model Q : 
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One may expect, that the projected model (Q) is mean- 
ingful only for the description of the low energy part of 
the eigenvalue spectrum - at higher energy the extended 
Hubbard rrt odpl may sh . 0W | additional structure due to the 
excitonk.se: idJto'feiohS, 'w'hich trte-hdithatecounted for in the 
projected rioSeT. 




FIG. 5. Lowest eigenvalues for all allowed momenta of 
12-site rings of extended Hubbard model ((7=10, V=5) and 
extended t — J model (J/t = 0.8) for 11 and 10 electrons (i.e. 
one and two holes in the half-filled band) . For the single hole, 
squares (triangles) denote states with S = 1/2 (S = 3/2), for 
two holes states with S = (S = 1). States of like spin have 
been grouped into 'bands' only to guide the eye. 

Figure || therefore compares the energies of the lowest 
eigenvalues of the extended Hubbard model with U/V=2, 
and the strong coupling model with J/t = 0.8 (we have 
used the spin correlation function at half-filling to 'gauge' 
the value of J; moreover, a nearest neighbor density re- 
pulsion with V = 5 has been added to the strong cou- 
pling model) . While the agreement between the two level 
schemes is not really pefect, there is an obvious corre- 
spondence. Deviations are largest for the 2 hole states 
with S=0 near 7r/3 = which probably indicates a 

somewhat higher 'holon' velocity in the actual Hubbard 
model. This may be related to the fact that U/V = 2 
is at the very borderline of stability against formation of 



a CDW ||, and the empty-double correlation function 
showed already some signature of the CDW. Apart from 
that, the strong coupling model reproduces the overall 
energy scales of the Hubbard model quite well, which 
indicates that for low energy excitations, despite its ap- 
parently unphysically large exchange constant, the strong 
coupling model may indeed provide the appropriate 'ef- 
fective' model. 

In summary, we have investigated the influence of inter- 
site Coulomb repulsion in the Hubbard model. For sin- 
gle band Hubbard models, our results suggest an en- 
hancement of the exchange constant, roughly as J(V) — 
4t 2 /(U-V), whereas the spin dynamics still remains con- 
sistent with that of an ideal Heisenberg antiferromagnet. 
The remarkable feature of the results is that this rescal- 
ing seems to hold up to the CDW instability, which opens 
the way to unexpectedly large values of the ratio J/t. In 
particular, for the strong correlation case close to a CDW 
instability in ID, our results suggest values of J compa- 
rable to the original hopping integral. Therefore, real 
materials which may be described by a strong correla- 
tion single band Hubbard model close to a CDW insta- 
bility may provide realizations of t— J-like models with 
'unphysically' large values of J/t w 1. For higher spatial 
dimensions, the critical value of U/V for the formation of 
the CDW state is smaller than in ID, which seems to pre- 
clude large J values due to inter-site repulsion. However, 
an alternative possibility here could be the suppression of 
the CDW due to a non-bipartite lattice structure, which 
can rule out the CDW even for large values of V; an inter- 
esting example here could be solid Cqq , which crystallizes 
in the fee structure. Auger spectroscopy fl gives values 
of U f» 1.5eV, whereas analysis of the exciton dispersion 
H in undoped Cqo suggests values of U — V as small as 
0.35eV, which is also consistent with other experimen- 
tal estimates Q. This is to be compared to the average 
bandwidth of the LUMO-derived bands of « 0.3eV M. 
While the inter-site Coulomb repulsion may be more ef- 
ficiently screened in the doped material than the on-site 
U, leading to an increase of U — V, one still may expect 
relatively large values of J. 
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